Neel transition, spin fluctuations, and pseudogap in underdoped cuprates by a 
Lorentz invariant four-fermion model in 2+1 dimensions 
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We show that the Neel transition and spin fluctuations near the Neel transition in planar cuprates 
can be described by an SU(2) invariant relativistic four-fermion model in 2+1 dimensions. Features 
of the pseudogap phenomenon are naturally described by the appearance of an anomalous dimension 
for the spinon propagator. 



The two dimensional one-band repulsive Hubbard 
model is considered as one of the best candidates 
for describing microscopically the planar cuprate high-Tc 
superconductors. The strong coupling limit of the Hub- 
bard model is equivalent to the t-J model In the t-J 
model at half-filling, on each site of a square lattice, on 
average a single electron interacts antiferromagnetically 
with its nearest neighbors and the system is an antiferro- 
magnetic (AF) insulator (a Mott insulator) described by 
Neel ordering. By introducing doped holes, thus charge 
carriers, on this lattice, the AF interaction is frustrated 
and a transition from the Neel ordered to the disordered, 
so-called spin-gapped phase or normal state occurs. For 
the Hubbard-Heisenberg model at and near half-filling 
Affleck and Marston showed, using a leading-order 
expansion, that the ground state is the tt-Aux phase for 
appropriate values of the hopping amplitude t, doping 
5, and AF interaction J ||^. The number N is the gen- 
eralization of the physical up-down spins, iV = 2, to iV 
types. In the tt-Aux phase the spinon spectrum has the 
dispersion 



2|x|a-\/cos^ aki + cos^ ak2 



where |x| is the absolute value of the 7r-fiux phase order 
parameter. This spectrum is gapless at the two Fermi 
vectors fp = (7r/2a, ±7r/2a) in the reduced Brillioun zone 
of the even and odd lattices with lattice spacing a. The 
linearization around these Fermi points gives a contin- 
uum (2+l)-dimensional massless Dirac theory describ- 
ing N flavors of four component Dirac spinors having a 
global U(2iV) symmetry [|, |]. At half-fiUing, the Dirac 
spectrum is isotropic and the flux-phase order parameter 
is equal to the so-called d-wave pairing order parameter 

|A| = |xl i- 

Recently Kim and Lee addressed the question, how the 
spin-gapped phase is connected to the Neel ordered phase 
at zero doping ||^. In their work the mean-field 7r-flux 
phase of Affleck and Marston is taken as the reference 
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state for describing the spin fluctuations around the AF 
Fermi points. By introducing gauge field fluctuations, 
enhancing AF correlations around the tt-Aux phase solu- 
tion, Kim and Lee propose, along the lines of Ref. that 
Neel ordering is described by dynamical symmetry break- 
ing (DSB) and mass generation in QED3. Neel ordering 
corresponds to the dynamically broken phase, which is 
characterized by a "mass gap" for the spinon spectrum 
and Nambu-Goldstone bosons as bound states of spinous 
and antispinons. These Nambu-Goldstone bosons are the 
massless AF spin waves. The disordered spin-gapped 
phase is equivalent to the subcritical symmetric phase 
and is characterized by the existence of massless spinous 
and unstable bound states as broad resonances. These 
resonances supposedly correspond to the spin excitations 
observed in the normal state and superconducting state 
of underdoped and optimally doped cuprates 

Although the physical picture sketched by Kim and 
Lee is plausible, it was pointed out that QED3 is not an 
appropriate model for describing unstable bound states 

. The main problem being that DSB in QED3 is not a 
phase transition of the second-order type, but a so-called 
conformal phase transition, which does not allow light 
unstable bound states in the symmetric phase [0. As 
an alternative for the gauge interactions, we propose that 
relevant, Lorentz- invariant four-fermion or four- Fermi 
(4F) interactions with an ultraviolet stable fixed point for 
the four-fermion coupling drive the AF ordering. The AF 
lattice Heisenberg interaction H = JX]<a; y> Sx ■ Sy, ex- 
panded around the two Fermi points, gives rise to SU(2) 
invariant attractive 4F terms in the action of the model. 
At sufficient strong AF coupling DSB occurs, giving rise 
to the Neel state. Despite the fact that the real temper- 
ature is not necessarily zero, the time dependent quan- 
tum fluctuations and ordering are given by a zero tem- 
perature 4F model. The mean-field tt-Aux phase local 
order parameter x describes the thermodynamic equi- 
librium state and is therefore time independent. In ad- 
dition, there is no need for a chemical potential in the 
proposed model, since the (nearly) half-Alling constraint 
has already been taken into account via the mean-Aeld 
equilibrium real bosonic Lagrange Aeld The present 
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idea is partly inspired by Ref. [Q, where it was suggested, 
by analyzing various lattice 4F operators, that the only 
relevant operators are those which are Lorentz invariant 
in the continuum. 

We adopt the spin liquid ansatz of Refs. |^; the spin 
liquid is described by the mean-field large- iV tt-Hux phase 
for low doping. The flux-phase order parameter \x\ de- 
pends on temperature and doping. The AF Heisenberg 
interaction is reinstated for this spin liquid. The fluctu- 
ations of the holes are ignored, and their effect is only 
included via their mean-field effect on reducing the AF 
exchange J to Jg = J(l — 5)^ Therefore, on a lattice 
with spacing a' — a/\/2, we consider the action 



dt 



.<x,y> 



with 



X! (ii^^t)i^9t - Xyx)ca{y,t) - Hi 



Hi 



(1) 



(2) 



<x,y> 



and where < x,y > denotes nearest neighbors on an 
isotropic cubic lattice. The index a i labels the spin 
components and the spin operator is 

where a are the Pauli matrices and c, are the spinon, 
antispinon operators. A particular representation for this 
Hermitian 7r-flux-phase hopping parameter Xyx is 

Xx±ai.x = i\x\i Xx±a2.x = MX\, 

with the nearest neighbor vectors ai — (a',0) and 02 = 
(0,a'). The low-energy behavior of the kinetic term in 
Eq. (|^) is known to be equivalent to a two-flavor massless 
Dirac theory with the action 



Sk 



dt J d^k{p4ida° + c{kij^ + k2j^)]i;a., (3) 

fe<A 



where c — 2\x\a is the "speed of light" and ip = ^'^7*^- 
The momentum cutoff A is naturally related to the lat- 
tice spacing via A ~ 7r/2a. The fields -0^ are four- 
component spinors, 

f Cela \ 

Co2a 
V Ce2a J 

where 1,2 labels the Fermi point and e, o labels fields 
on the even and odd lattices, respectively. The 4x47 
matrices satisfy a Clifford algebra corresponding to the 
Minkowskian metric g^"^ = diag(l, — 1, — 1). The follow- 
ing representation for the 7 matrices has been chosen: 



V'i — {cli^ C^i^ c],2q Cg2Q^ , (4) 



7 







1 iai 






1' = 


\ 







ia2 

-i(T2 



where ai are Pauli matrices, acting on the even and odd 
site fermion operators. 

Expanding the AF Heisenberg interaction (^) around 
the two Fermi points, we obtain 



Hi = 4Jea" j (2^)"<5(pi + fci - - P2) 

Pl,P2,fcl,fc2<A 

[4la(Pl'*)Cel/3(P2,t) +42a(Pl,^)Ce2/3(P2,i)] 
X [4l7(^l'*)Col5(fc2,t) + cl2^{ki,t)Co2s{k2,t)] 
X (SasSj/l - Sa/sSysf^), (5) 

where the SU(2) Fierz identity has been used, cfafs ■ a-ys — 
2SasS-yf3 — SapS-fp. Using Eq. (jj), Eq. (||) can be written 
as 

Hi = ~Jea^ J {2TTf5{pi+ki-k2-p2) 

pi,P2,fel,fe2<A 



(6) 

Subsequently, it is straightforward to show that Eq. (||) 
together with Eq. gives rise to the action S = Sk + Si, 
with 



Si = .W J dt J {2nY5{pi +ki-k2- P2) 

Pi,P2,ki,k2<A 

(7) 

with the generators of the SU(2) symmetry with 
Tr (r^r^) = 6-^^ (r^ = a^/V2), and V^r^^V = 
'^a f3^o''''a0'^l3- ^hc action (|^) is invariant under global 
SU(2)xU(l) corresponding to the spin orientation sym- 
metry and total spin conservation. Moreover the action is 
invariant under the discrete transformations; space reflec- 
tion, parity, and the combined CT (charge-conjugation 
and time-reversal) transformations. Naturally, it is in- 
variant under continuous rotations in space. However, 
the terms of the form -0^?/' are not invariant under Lorentz 
boosts and therefore the action is not relativistic invari- 
ant. 

In what follows, we shall show that the action (|^) 
lies in the same universality class as that of a Lorentz- 
covariant SU(2) invariant (2+l)-dimensional 4F model 
for two massless fermion flavors: 



/ 



d-^x 



cG 



(8) 



where xq — ct, d ^ j^d/dx'^, and where G is an at- 
tractive four-fermion coupling, cG/2 ~ J^a^. The uni- 
versality only holds close to a critical point or ultraviolet 
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stable fixed point Gc of a second-order phase transition of 
Eq. (^) at which the SU(2) symmetry is dynamically bro- 
ken to a U(l) symmetry. There is a dynamical generation 
of a parity-conserving "mass term" nig connected with 
the appearance of a nonzero vacuum condensate {-ipT^ip). 
This Lorentz-invariant condensate describes a staggered 
spin expectation value giving rise to AF ordering ^, 
In the broken phase, two massless pions or spin waves 
appear as Goldstone bosons 0. 

The spinon propagator Sap can be written as 



'q/3 



{p)p^ — ml 



with Minskowskian momentum p^ ^ Pq ^ P' P and where 
A(p) is the fermion wave function. In the Hartree-Fock 
approximation, the equation for gets contribution 
only from the tad pole diagram and the fermion wave 
function A{p) ~ 1. Setting c ~ 1, the gap equation for 
TO,, reads 



Gi 



(Pp Arris 



M 



(27r)3 p2 



with the subscript M denoting the Minskowskian metric 
with cutoff A. This gap equation gives rise to the familiar 
critical coupling gi^'^ = 2GcA/7r^ = 1. Above the critical 
couphng g — 2GA/tt^ > ^c^^ the SU(2) spin symmetry is 
broken and a {-ipT^tp) condensate is formed. 

Now let us show that the nonrelativistic interaction 
terms in Eq. (Q) are irrelevant close to gi^'' ■ We in- 
vestigate the generation of a mass to„ connected with 
the uniform spin expectation value (ip^T^ip) @- Iii the 
Hartrcc-Fock approximation, keeping only the Lorcntz- 
noninvariant terms, the spinon propagator is of the form 



Saflip) 



(pSafj " niu-ioTafs) 

{[po~i-ir+'mjV2f-pl 



with a — \,2 {a =t,i)- The gap equation for to„ then 
reads 



4Gi 



M 



d^p 

(27r)3 



Naturally, the Hartree-Fock approximation ignores 4F 
fluctuations and therefore gives a rather crude descrip- 
tion of the critical behavior. Nevertheless the irrelevance 
of the non-Lorentz-covariant terms in Eq. ^ to gi^^ can 
be demonstrated in more advanced approximations, such 
as the 1/A'' expansion. In particular, the SU(2) Heisen- 
berg antiferromagnets can be generalized to SU(7V) and 
studied in the 1/N expansion This is analogous to 
the SU(iV) generalization of Eq. (^). Contrary to models 
of the Gross-Neveu type , such an expansion resem- 
bles the topological 1/N expansion of 't Hooft |l3j, cor- 
responding to an expansion in planar Feynman diagrams 
instead of Fermion loops. 

Recently, the generalization of Eq. (||) to N fermion 
flavors with a SU(A^) xU(l) invariant 4F potential con- 
taining iV^ — 1 terms has been studied in the planar large- 
N approximation After a Hubbard-Stratonovich 

transformation (a"^ = —G'ipT^ip), the action (H) can be 
expressed as 



S 



A A 



(9) 



In the SU(A^) case, the spin label A runs from 1 to N"^ — 
1 and the fermion label from 1 to N . Thus there are 
N"^ — 1 composite operators, giving rise to N"^ — 1 spin 
propagators (connected), iA^{q) = {aA{q)crA{—q))c, and 
iV spinon propagators, iiS'(p) = {ipip)ip{—p))- This allows 
for a 't Hooft topological 1/N expansion as is conjectured 
in Rcf. Moreover, it was argued in Ref. |Q that 

the leading large-TV or planar approximation reduces to 
the ladder approximation for the so-called Yukawa vertex 
{i.e., r^{k,p) = T^^l). The Yukawa vertex is the fully 
amputated three-point interaction vertex for the action 
(^). In the ladder approximation, the Schwinger-Dyson 
equations for the spin propagators and the spinon 
propagators S form a closed set, as depicted in Fig. ^ 



-O- 




The bifurcation equation is 
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M (2^)^ (pI 
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PI 



2^2 ' 



pI) 



giving rise to a critical coupling g, 



(u) _ 



2GcA/7r2 ^ 3. 



This shows that when g is close to the critical value g, 



it is far (\g — g^c^\/ g — 2) from the critical regime of g*-"-*. 
Upon increasing g it first encounters the critical value 
for staggered magnetization. Hence, the nonrelativistic 
terms in Eq. are irrelevant in that area of the coupling 
constant space. 



FIG. 1: The leading large-A'^ truncation or ladder approxi- 
mation for the Schwinger-Dyson equations for the spin- wave 
propagator Act(p) (dashed line with blob) and the spinon 
propagator S'(p) (solid line with blob). 

For the physically relevant case [N — 2), the prop- 
agators A;^, A^ of the auxiliary fields cti, describe 
the Goldstone modes, whereas the field (T3 acquires 
a nonzero, but Lorentz-covariant, vacuum expectation 
value in the broken phase ((era) 7^ 0, g > g^c^^, describing 
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the Neel state, corresponding to the symmetry-breaking 
pattern SU(2)^U(1). In the subcritical region (g < gi^^) 
the propagators of all three auxiliary fields describe Gold- 
stone precursor modes that come down in energy as the 
transition is approached. Hence, the staggered spin fluc- 
tuations or spin waves are described by the propagators 
of the auxiliary fields a a , which become light close to the 
critical point g = ^ic"^ . 

In Ref. ||lj], the Schwinger-Dyson equations repre- 
sented in Fig. 1^ were solved. It was shown that the 
fermion or spinon propagator S acquires an anoma- 
lous dimension C via the fermion wave function A{p) ^ 
(A/p)^, so that at the critical coupling G = Gc the 
fermion (spinon) propagator scales as 



(10) 



The dependence of the anomalous dimension ( on N is 
determined, and for iV = 2, we have C. ~ 0.21 [Q. More- 
over the model described in Ref. Jl^ , with the dimension- 
less 4F coupling g = 2GA/tt'^, has an ultraviolet stable 
fixed point ai g — gc — 1 + 2^. The appearance of a pos- 
itive anomalous dimension for S is considered to provide 

0- 



a description of the pseudogap phenomenon llq, 16 



An important point is whether the critical coupling 
1 -|- 2C is in agreement with the estimations for the 



JeCl 



2|x|a, A: 



physical parameters. Since Gc/2 
7r/2a, we have that g ~ Ji,/n\x\- With the estimation 
f/c"'' — 1.4, we obtain that a value \x\ ~ 0.23Je would 
get us close to criticality. This is remarkably close to the 
mean-field value |x| ~ 0.24J given in Ref. Q. 

In the ladder approximation (Fig. |l|), the connected 
propagator of the cta field in momentum space reads 



A 



A-\q) = -±+t I d^kTY [T^S{k + q)T^S{k)], 



1 



G 



M 



where S{k) is the full spinon propagator, given by 
Eq. (|l^), see also Fig. In the subcritical or SU(iV) 
symmetric regime, we can take S{k) = k/[k'^A{k)]. The 
integral can be performed, and A^(g) has the following 
scaling form for j^j ^ A (g^ = Qa ^ Qi ^ li)- 



A^(g) 



C (-AV92)C+l/2 



A 



1 + (-m2/q2)C+i/2 



(11) 



where C is some flavor dependent positive constant p4[ . 
The mass Too- (spin-wave stiffness) denotes the position 
of the resonance peak given by the imaginary part of 
Eq. ( pj]), and plays the role of the inverse correlation 
lengthl. 



1/(1+20 



(12) 



From these expressions, it follows that the critical expo- 
nents 77, v, and 7 are 77 = 1 — 2^, 1/ = 1/(1-1- 2^), and 



7 = 1. These exponents satisfy the three-dimensional hy- 
perscaling equations jl^ . Moreover, in Lorentz-invariant 
field theories, the scaling of the energy equals the scaling 
of momentum, and consequently the dynamical scaling 
exponent z = 1. 

In experiments, the isotropic dynamical susceptibility 
x"(g) is measured with q = {u>,q). In Ref. |^ mag- 
netic resonances were observed in underdoped and op- 
timally doped YBa2Cu306+K, with the famous 41-meV 
peak at optimal doping. For lower doping the resonance 
peak shifts to lower energies. It is tempting to assume 
that, along the fines of Refs. ^, these resonances might 
be described by the Dirac models, see Fig. ^. However, 




FIG. 2: The spin susceptibility or spin-correlation function 
Im[ACT(p)] vs p/rria, with the peak renormalized at unity, 
C = 0.21 [ri = 0.58). 

optimal doping {6 ~ 0.2) is rather far from Neel dop- 
ing, which is one order of magnitude less, 5c — 0.02 [p^ . 
Therefore the question is whether we are "sufficiently" 
close to the scaling region of the Neel transition. 

Let us compare a couple of (2-|-l)-dimensional Lorentz- 
covariant quantum field models, capable of describing AF 
fluctuation and Neel transition. For instance, Kwon jl^ 
presented a nodal d-wave spin liquid model for the Neel 
transition, which has the universality class of the Gross- 
Neveu model. This model has a single AF order parame- 
ter, corresponding to the singlet composite order param- 
eter {^ip), with a single two-component Dirac fermion 
{N =1). The universality class of the model presented 
in this paper {i.e., Eq. mf) deviates considerably from 
the Gross-Neveu model [Q. Although the anomalous 
dimension rj for both models is comparable: i] « 0.58 for 
the present model, and 77 = 16/{STr^N) « 0.54 {N = 1) 
for the Gross-Neveu model, the main difference is in the 
value of the anomalous dimension of the fermion propa- 
gator, which for the present model is C ~ 0.21, whereas 
for the Gross-Neveu model it is 77 = 2/{3tt'^N) « 0.07 
{N — 1). Another approach was adopted in Ref. [ pO[ , 
where the antiferromagnetic correlations and the Neel 
transition are described by the nonlinear sigma model 
in the large-TV expansion. In that work, the order pa- 
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rameter is a three-component one, corresponding to the 
three spin-1 components, and no Dirac fermions are taken 
into account. The anomalous dimension ij for the spin- 
waves turned out to be 77 = 8/(37r2iV) « 0.09 {N = 3), 
which is considerably smaller than 77 for the two above 
mentioned 4F Dirac models. In this paper, we have a 
three-component order parameter (the composite spin 
degrees) and two flavors {N = 2) of four-component 
Dirac fermions. To determine which universality class de- 
scribes the Neel transition, the low doping region 6 ^ 6c 
needs to be examined in more detail experimentally. 

Since in the subcritical phase, all three correlation 
functions of the staggered spin components are degen- 
erate, the imaginary part of the correlation function 
A;^(g) is directly proportional to the so-called odd acous- 
tic mode of x" (9) fl^- The spin-correlation function only 
gets low-energy contributions from the staggered spin op- 
erators; 

{Sa{p)Sb{-p)) (X {aA{p)<JB{~p))o,SABiA^{p). 

Moreover, since g is proportional to Je, g reduces when 
doping is increased The critical coupling Oe''^ of the 
Neel transition is a quantum critical point |21[| and cor- 
responds to a critical Neel doping (5c ^ 1- Equation (|l^ ) 
gives the relation between the position of the peak of the 
magnetic resonance and the doping rate 6. Assuming Sc 
is sufficiently close to zero, we obtain that the resonance 
peak is linearly proportional to doping (for small dop- 
ing rates with 6 > 6c)- Consequently, the peak position 
moves to lower energies when doping is reduced; near dc 
the peak height diverges. For doping values 5 < Sc spin 
waves appear. 

In summary, we have shown that a (2-|-l)-dimensional 
Lorentz-invariant 4F model with a global spin SU(2) sym- 
metry describes the low-energy time-dependent "quasi- 
particle" spin excitations of the t-J model near the AF 
wave vector. The spin excitations are given in terms 
of quantum fluctuations around the mean-field 7r-fiux 
phase. The Neel transition is described as the DSB of 
SU(2)— >U(1) in the model. The magnitude of the critical 
coupling gc of the 4F model turned out to be in good 
agreement with the input parameters J and which 
define the fiux-phase spin liquid. Nevertheless, the ques- 
tion of the precise effects of hole doping on the AF inter- 
action and the anisotropy of the Dirac spectrum is left 
open. For the future, it would be interesting to include 
the contribution of the slave bosons (holons) {e.g, see 
Refs. 1^, 18 ) in order to take into account anisotropy and 
to determine the effective AF coupling Jg. We demon- 
strated that the unstable spin modes found in experi- 



ments might be well described by the Nambu-Goldstone 
boson precursor modes in the subcritical region. Fur- 
thermore, the spinon propagator acquires an anomalous 
dimension, but remains gapless near the AF wave vec- 
tor in the normal state. The appearance of a spinon 
anomalous wave function gives a natural description of 
the pseudogap phenomenon. 

The author thanks V.P. Gusynin for useful suggestions. 
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